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Abstract 

We study the quark mass expansion of the axial-vector coupling constant gA of the 
nucleon. The aim is to explore the feasibility of chiral effective field theory methods 
for extrapolation of lattice QCD results — so far determined at relatively large quark 
masses corresponding to pion masses m n > 0.6 GeV — down to physical values of m^. 
We compare two versions of non-relativistic chiral effective field theory: One scheme 
restricted to pion and nucleon degrees of freedom only, and an alternative approach 
which incorporates explicit A(1230) resonance degrees of freedom. It turns out that, 
in order to approach the physical value of gA in a leading-one-loop calculation, the 
inclusion of the explicit A(1230) degrees of freedom is crucial. With information 
on important higher order couplings constrained from analyses of the ttN — > ttttN 
reaction, a chiral extrapolation function gAim-w) is obtained, which works well from 
the chiral limit across the physical point into the region of present lattice data. The 
resulting enhancement of gAim-w) near the physical pion mass is found to arise from 
an interplay between long- and short- distance physics. 
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1 Introduction 



Lattice QCD is developing into a powerful tool for studying the structure of nucleons [UI2]- 
In practice, however, these computations are so far limited to relatively large quark masses. 
The typical "light" quark masses that can be managed on the lattice are usually at least ten 
times larger than the u- and d- quark masses, m Uj d — 10 MeV, determined at renormalization 
scales around 1 GeV. The corresponding masses of a pion on the lattice — defined as the 
lowest lying 0~ state in the simulation — are larger than 0.5 GeV, quite far above the physical 
value m n ~ 0.14 GeV. 

In the limit m u ^ — > 0, on the other hand, low-energy 2-flavor QCD displays a sponta- 
neously broken chiral SU (2) x SU (2) symmetry. Pions as massless Goldstone bosons are 
the relevant degrees of freedom of the (effective) theory. Their coupling to baryons is also 
subject to the rules imposed by chiral symmetry. It produces the pion cloud of the nucleon, 
an important low-energy, long- wavelength aspect of nucleon structure, cf. [2j and references 
therein. 

For baryon properties the interpolation between lattice results obtained at relatively large 
pion masses and actual observables determined at the physical m n has become an issue of 
great interest recently. Chiral effective field theory (ChEFT) can, in principle, provide such 
extrapolations. First steps in this direction were made by Leinweber et al. jS] for the case of 
the magnetic moments. They used Pade approximants based on the leading dependence on 
the pion mass as dictated by chiral symmetry. A perturbative analysis of nucleon magnetic 
moments, using ChEFT with explicit inclusion of the A(1230) isobar, has been performed 
in jlj. It turned out to be of crucial importance to promote the isovector couplings of the A 
to leading order in the power counting characteristic of chiral effective theories. Truncated 
versions of chiral perturbation theory usually "freeze" the A and relegate its effects to higher 
order counterterms, with the consequence that the range of convergence of this truncated 
theory is often quite limited. The importance of the A as an explicit degree of freedom in 
various aspects of nucleon substructure has long been known. It has played a prominent 
part in chiral models of the nucleon [21 E] and earlier studies based on current algebra [S] . 

In the present paper 2 we investigate the quark mass dependence of g^, the axial- vector 
coupling constant of the nucleon. This is a key quantity for our understanding of the 
nucleon's chiral structure [2]. Furthermore, lattice QCD determinations of qa are progressing 
to the point where such an investigation is in demand for interpolating between lattice data 
and the actual observable [TJ |H] . 

Of course plenty of ChEFT calculations pertaining to this quantity exist in the literature 
jOj, and the special importance of intermediate spin 3/2 resonance contributions for the axial 
properties of baryons has been recognized a long time ago. Recall, for example, the Adler- 
Weisberger sum rule [TU] : 

2 f 2 r°° du f m 2 \ 

A = ! + r~ 2 f K+» ~ VpM] + O ^f- ■ (1) 

It relates the surplus of the axial- vector coupling constant beyond its "trivial" value qa — 1 
2 Some aspects of this work have already been reported in [Jj. 
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(for a point-like, structureless nucleon) to the excess of the n + p cross section over the it p 
cross section, a feature that is dominated by A (1230) resonance excitation. 

In ChEFT qa is calculated from the response of a baryon to the presence of an ex- 
ternal background field with axial quantum number. In contrast to many earlier ChEFT 
calculations — e.g. the pioneering ones of refs. [TT| IT2"] — we are not only interested in the non- 
analytic quark mass behaviour generated by the Goldstone boson cloud around the baryon. 
In order to analyze the quark mass dependence of a baryon observable in a quantitative way 
we follow the same philosophy as spelled out in ref . [3] . We specify a power counting scheme 
plus a certain order in that scheme and then systematically evaluate all contributions — 
short- and long-distance — to that order. It turns out that the role of short- distance physics 
parameterized as higher order operators in chiral effective field theory is crucial for an un- 
derstanding of the chiral extrapolation function ^(m,,-). We utilize the so-called SU(2) 
"small scale expansion" of ref. ^2] to leading-one-loop order (0(e 3 )), which includes explicit 
nucleon and A(1230) degrees of freedom — some details regarding the formalism are summa- 
rized in sectional In fact, all the long- and short-distance contributions considered in our 
analysis have already been found (as a by-product) in ref.^3], where among other topics the 
(small) momentum dependence of the axial form factors of the nucleon was studied. Here we 
reconsider these results and focus on the quark mass dependence of g^- Unlike the situation 
in the vector current sector, where a modified power counting had to be developed jl] so 
that important effects of the A could be captured already in leading-one- loop order, this 
turns out not to be necessary for axial current effects related to the A(1230). For the axial 
structure of baryons — which is the topic of the present analysis — the standard counting of 
ChEFT, as employed in refs. |T3*1 114j. is sufficient to obtain A(1230) induced quark mass 
dependence already at the leading-one-loop level. Throughout this work we can therefore 
apply standard ( "naive" ) counting rules. 

Historically, attempts to obtain a chiral extrapolation function for based on the known 
leading-non-analytic (LNA) quark mass dependence in combination with a phenomenolog- 
ical (quark mass dependent) regularization procedure did not yield satisfactory results, 
displaying axial couplings less than unity at the physical point [T5] . In our opinion one 
should not have been surprised about such a failure, as the LNA quark mass term presum- 
ably is only dominant for quark masses near the extreme chiral limit of the theory. Such 
a feature was indeed observed in the analysis of the anomalous magnetic moment of the 
nucleon jl]. Recently, Detmold et al. ^H] in their analysis of moments of polarized DIS 
structure functions, found an improved extrapolation formula for (yu^vr) utilizing a chiral 
quark model which also allows for contributions from intermediate A(1230) states How- 
ever, the resulting extrapolation function — which has most of the A(1230) related couplings 
fixed from SU(6) symmetry — still does not provide for an enhancement of ^(m„-) near the 
physical point, which would at the same time connect with the lattice results. In this work 
we go beyond the existing analyses by performing a systematic perturbative calculation of 
all short- and long-distance contributions which are allowed to leading-one-loop order in 
chiral power counting. 

This paper is organized as follows. In the next section we give a brief summary of the 
ChEFT formalism utilized here. In section El we present the leading-one-loop analytic result 
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for both in a scheme without and in one with explicit A(1230) degrees of freedom. In 
section |3] we compare these analytic results with recent data from the QCDSF collaboration 
[TT] . employing three different fit procedures. Finally, in section El we check the stability of 
the obtained results by including the physically known information about qa as an additional 
constraint in the fit. The resulting chiral extrapolation function ^(^tt) provides a sensible 
description of the quark mass dependence from the chiral limit over the physical point out 
to the present lattice data. We conclude with an outlook on further investigations required 
for a better understanding of gA- 



2 Chiral Effective Field Theory Input 

Our analysis of gA is based on chiral effective field theory with two light flavors (u, d). We 
work to leading-one- loop order. All other quark degrees of freedom are integrated out, 
leaving their marks only in slightly shifted values of couplings sensitive to short distance 
dynamics. This effective field theory of low-energy QCD acts with pions as Goldstone bosons 
of the spontaneously broken chiral SU(2) x SU(2) symmetry of QCD. In addition the chiral 
symmetry of the QCD Lagrangian with two light flavours is explicitly broken by the small 
u— and d— quark masses which are treated as a perturbation. They shift the massless 
Goldstone bosons to the physical pions with mass m n and lead to a string of quark mass 
dependent operators which turn out to be of crucial importance for the understanding of 
chiral extrapolation functions. At present we ignore effects of isospin breaking and work 
with degenerate masses, rh = (m u + m < i)/2, for the up and down quarks. 

We utilize two versions of ChEFT for our analysis of gA- The first scheme is SU(2) 
heavy-baryon chiral perturbation theory (HBChPT), which involves only the Goldstone 
boson modes (pions) and spin 1/2 matter fields (nucleons), taken up to order p 3 . The 
second scheme is the SU(2) small scale expansion (SSE) which in addition also involves 
explicit spin 3/2 matter fields (four A(1230) states) and their interactions with Goldstone 
bosons, taken up to order e 3 in the "small scale" e which now includes also the non-zero 
A-iV mass difference |18j . 

For completeness — and for the proper definition of our couplings — we give the effective 
Lagrangian required for this leading-one-loop analysis of g^- 

The well-known leading order pion Lagrangian C^l can be found in ^2] • The leading order 
tiN, iiNA and nA Lagrangians 

£(J) = N v [ivD + g° A S-u]N v , 
£$ A = ffc A wiN v + h.c. , 

£ W = Tt[zt 3/2 A -tvDV- gi S-u5^g, u T» (3) 

are specified in detail in • Here we only note that N v (T/ 4 ) represents the non-relativistic 
spin-1/2 nucleon (spin 3/2 delta) field, with £}"L 3 / 2 denoting the isospin 3/2 projection op- 
erator. The chiral tensors m m , encode couplings to (external) axial sources with NN-, 
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iVA- and A A- axial coupling constants g A , ca and gi, denned in the chiral limit. Fur- 
thermore, D^, D l j denote the chiral covariant derivatives of the nucleon, respectively the 



Delta. A corresponds to the iVA mass splitting in the chiral limit. The two four- vectors 
v^i, Sfj, — related to kinematics and spin — are discussed in the Appendix. It can be clearly 
seen from Eq. (jH)) that all couplings considered as leading order follow the standard counting 
ofChEFT. 

While the leading order Lagrangians given above are well-known in the literature we also 
want to discuss the less-known higher order couplings required for a complete leading-one- 
loop calculation of gA both in HBChPT and in SSE. We utilize 
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Here Mo denotes the nucleon mass in the chiral limit, whereas the chiral tensor \+ encodes 
quark-mass dependent short distance physics. For the higher order couplings Bi we follow 
the nomenclature of ref. ^3] where a complete set of counter terms is listed, appropriate for 
renormalization of both leading one-loop HBChPT and SSE calculations. We note that B 3Q 
and B 3 i are identically zero in HBChPT, but are required for leading-one-loop calculations 
in SSE. In general all Bi have a finite regularization scale A dependent part -B[(A) as well 
as an infinite part: 



Bi = Bl{\) 



0i 



(4tt/ w ) 



16tt 2 L 



(5) 



In the infinite part denotes the HBChPT (respectively SSE) beta-function associated 
with this counterterm, whereas the infinities are encoded in the quantity L, discussed in the 
Appendix. f w is the pion-decay constant. 



3 Analytic Results 



We now proceed to a discussion of the results for the axial coupling of the nucleon expressed 
as a function of the pion mass gA{ m n), first for the HBChPT and then for the SSE case. The 
leading one-loop HBChPT result for the quark mass dependence of gA is known for a long 
time |HJ. Evaluating diagrams 1-5 of FigGJand projecting out the axial coupling constant 
gives 



g A (ml) = g A Z N + AmlBl{\) 



m; 



32ttV= 



ln^ + (^) 3 



0{p A ) , (6) 
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Figure 1: Diagrams contributing to the nucleon axial- vector coupling constant g A at leading- 
one-loop order. The wiggly line denotes an external (isovector) axial- vector background field, 
interacting with a nucleon (solid lines). 



where g° A again denotes the axial- vector coupling constant in the chiral limit. The Z-factor 
of the nucleon to leading-one-loop order reads 



3(g° A ) 2 ml ( m Tl 



Z N (mi) = 1 - ™ 2f2 * (1 + 3 In ) - 8m^ (A) + 0(p*) . 



Following ref . [20] . we note that this Z-factor is finite since we have not transformed away 
the counter term B 2 o- It is well-known, however, that the finite part of a low-energy counter 
term such as B 2 q cannot be observed independently (of -Bg(A) in our case). To order p 3 
in HBChPT, g A {m\) therefore depends on only two unknown parameters, g\ and C[X) = 
Bg(X) — 2g A B'! 20 {\) ) for each choice of the regularization scale A: 

9 H A B {ml) = 9 A -^rnl + i{c HB (X) +1 HB \n^}rnl + 0^), (8) 



with 



1 HB 



(9°a) 3 + 19°a 



(9) 



We expect Eq.flBJ to hold for sufficiently small quark masses, close to the chiral limit. All 
observable short distance dynamics is collected in C(A). The procedure outlined here is 
well-defined in the sense that the scale-dependent logarithm in Eq.flBJ) cooperates with the 
scale-dependent coefficient C(A) in just such a way that the overall sum is scale-independent. 

Next we turn to the result of the SSE approach which includes explicit A (1230) degrees 
of freedom. The bare result depending on all four counterterms of Eq. (j3J) has been given in 
ref. [2] . It originates from diagrams 1-8 of Figd We note that diagrams 2 and 3 in principle 
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also exist when the intermediate nucleon is replaced by a spin 3/2 baryon. However, the 
leading order NA transition Lagrangian of Eq.© does not include operators connecting 
iVA via an even number of axial fields. The corresponding diagrams are therefore of higher 
order than the ones considered here at the leading-one-loop level. Extending the work of 
ref.jH] we also impose the constraint of decoupling which implies that, in the chiral limit, 
effects of the A(1230) can be absorbed in contact interactions. This makes sure that the 
addition of explicit Delta degrees of freedom maintains the correct m, — > limit of g A - 
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and the leading-one-loop SSE Z-factor given in the Appendix we find the chiral limit be- 
haviour 
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We emphasize the point that, while the result is unique, the separation into B^ Q and 
B^ can be done in several possible ways. Comparing with Eq.fjHJ) we note that with our 
parameter choice of Eq.ljlUJl we have restored g° A as the chiral limit value of gAim 2 )- Further- 
more, we note that the leading quark mass dependence starting proportional to m\ \nm n is 
not modified by the addition of explicit A(1230) degrees of freedom, again consistent with 
chiral symmetry. 

Having taken care of the decoupling constraints we now present the complete leading- 
one-loop result for g A in the Small Scale Expansion: 
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Table 1: Values of the parameters which are taken at their physical value 



Parameter 


Physical Value 


fir * fir 


0.0924 GeV 


C° A -> C A 


1.125 


A -> Re[A] 


0.2711 GeV 



We note that the first line in Eq. (jl2j) has the same structure as the heavy baryon result 
of Eq.(jHJ). New structures appearing in the leading-one- loop SSE result are the terms pro- 
portional to m 3 as well as the logarithms depending explicitly on the N-A mass splitting 
value. For the chiral extrapolation of g\ we will utilize the full analytic form as given in 
Eq. (J12|) . The much simpler chiral limit form Eq. (jllj) will only be used to constrain coupling 
constants. 

4 Chiral Extrapolation to Lattice Data 
4.1 General Remarks 

We now turn to a numerical evaluation of the two analytic extrapolation formulae for qa 
considered here. In principle all couplings and masses — aside from m n — on the left hand 
side of Eqs. (j8ll2|) are understood to denote their values in the chiral limit. However, for 
many of them the chiral limit values are only poorly known. For the pion decay constant 
we use its physical value because the difference between this and /° is known to be only 
a few percent. For ca and Ao — following the reasoning laid out in ref.[l] — we also use the 
empirical values, specified in table ^ 

Choosing A = 1 GeV without loss of generality in leading-one- loop SSE (HBChPT), we 
are then left with three (two) poorly known couplings: g% C SSE (1 GeV), g±, or g% C HB (1 GeV), 
respectively. To determine these short- distance physics parameters of interest, we utilize 
the (quenched) simulation data of the QCDSF collaboration ^7] in the mass region 
0.3 GeV 2 < m 2 < 0.6 GeV 2 , containing 5 data points. Hereby we are working under the 
assumption — following again the reasoning of ref. jlj — that for lattice pion masses larger 
than 600 MeV the effects of "quenching" can be neglected. This intuitive expectation has 
recently been put on firm ground by the analysis of ref. jH] . It was concluded that for present 
state-of-art simulation data for moments of nucleon structure functions (such as g^) indeed 
no significant differences between quenched and fully dynamical QCD simulations could be 
found for > 500 MeV. However, the difference 3 between "quenched" and fully-dynamical 
QCD simulations should become visible at lower pion masses. 

While we are confident that (un) quenching effects are small in the mass range considered 

3 In quenched QCD the axial coupling of the nucleon develops a chiral singularity gAijn^ — > 0) ~ logm^ 
in the chiral limit m contrast to the finite value <?a(tott — > 0) — ► g\ in full QCD. 
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here possible corrections arising from the finite simulation volume could turn out to be 
important for the exact position of the lattice data. The five data points we are utilizing 
here actually result from three different lattice spacings which, however, roughly correspond 
to the same physical volume 4 of (L=1.45. . . 1.60 fm) 3 Given this data situation no 

conclusions on finite volume dependence can be drawn from the QCDSF data. In ref.[2HJ for 
a simulation utilizing domain wall fermions it has been reported that qa is a quantity which 
can be sensitive to such corrections. We note that the data analyzed here are consistent with 
one another (c.f. Fig|2J) and are basically flat in their m n dependence. The data shown in 
ref . are also rather flat for the largest volume considered there, albeit the overall location 
of the plateau seems to lie slightly higher than the QCDSF data discussed here. At present 
no systematic study of these effects exists and we therefore assume that possible shifts in 
the data due to the finite simulation volume fall within the size of the error bars of the data. 
Once a larger amount of data for qa^-k) are obtained at different lattice sizes, we will take 
up this issue again to recheck our analysis. 

For the purpose of our numerical analysis the qa simulation data thus only provide the 
constraint of specifying the location of the typical "data plateau" (e.g. see FigJ2J), which is 
quite analogous to the situation in the case of the anomalous magnetic moments [3]. Once 
more, we emphasize the empirical observation that any sensible chiral extrapolation formula 
must have enough structure to reproduce such a plateau, at least over a certain range in m^. 
In essence, this implies a strong constraint: different terms in the quark mass expansion of 
qa must co-operate in such a way that the pion cloud effects are basically counterbalanced 
by short-distance physics once the pion mass exceeds 0.6 GeV. Given that the available 
lattice data are still restricted to relatively large quark/pion masses, can one nevertheless 
utilize this information to determine the chiral extrapolation function? In order to answer 
this question we perform three different fit procedures of varying sophistication. 

4.2 Free Fits 

In the first round of fits we want to explore whether the available lattice data are sufficient 
to constrain the unknown parameters g A , C SSE (X), gi for SSE, or respectively, g A , C HB (\) 
for HBChPT. The numerical values obtained are given in table 121 with the labels fit la, 
fit lb. The results are shown in FigJ2] via the short-dashed (HBChPT) and long-dashed 
(SSE) curves. While the \ 2 P er degree of freedom (c.f. table El) for either curve is small, 
clearly both extrapolation functions have to be considered unphysical, as indicated by the 
small values of g A violating the Adler-Weisberger sum rule as well as the physical constraint 
gPhys. _ ^2670(30) obtained from neutron beta decay [23]. 111 addition, the SSE curve 
has a value of the axial A A coupling constant g\ « 0, whereas we expect g\ > g A from 
SU(6) symmetry considerations. We note that the resulting (unphysical) HBChPT curve is 
quite similar to one found in ref.^2], where the chiral logarithm with r ) HB of Eq.(|H|) plus a 
phenomenological (quark mass dependent) cutoff procedure parameterizing short-distance 

4 Whilc this volume at first glance appears rather small to study the physics of an extended object like 
a nucleon, we note that the nucleon radii are considerably smaller at these large quark masses than at the 
physical point. See ref. |22j for an analysis of the quark mass dependence of the electromagnetic radii. 
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Table 2: Fit results discussed in the text 



Fit 


n° 
9a 


V I /\ — Xvjrc V J vjc V 


yi 




v 2 /H n f 

A /U.O.I. 


la 


U. i 1 ± U.U4 


n in in HQ 




Q 
O 


n Q7 
U.o / 


lb 


0.78 ±0.04 


±1.06 ±0.08 


0.0 ±0.1 


2 


0.55 


II 


0.94 ±0.04 


-0.25 ±0.04 


9/5 g° A 


3 


0.39 


III 


1.12. ..1.26 


-2.2... -4.6 


4.3... 6. 6 


3 


0.48... 0.60 


IVa 


1.36 ±0.07 


-5.3 ± 1.2 








IVb 


1.21 ±0.01 


-3.4 ±0.4 


5.6 ±0.5 


3 


0.54 


IVc 


1.21 


-5.3 ± 1.2 









physics was employed. We note that both the HBChPT as well as the SSE "free fits" 
prefer a small positive value for the short distance physics parameterized via C(X) — at a 
regularization scale of A = 1 GeV. The reason for this unsatisfactory situation is that the 
present lattice data are not yet sufficiently accurate to obtain a realistic value for C, as we 
will see in the next two sections: additional physics constraints need to be invoked in order 
to limit the range of C. 

4.3 SU(6) Constrained Fit 

One possible conclusion from the mismatch between the chiral extrapolation and present 
lattice data could be that the extrapolation formulae (jHI) and (112)) are just too simplistic 
to cover such a large region in (for example, this was the case with the corresponding 
extrapolation functions of schemes A or B of the anomalous magnetic moments in ref.jl]). 
In the following we argue that this is indeed true for the leading-one-loop HBChPT result 
of Eq.(jH|), whereas the corresponding SSE extrapolation function Eq. (jl2)l . with the A(1230) 
added as an explicit degree of freedom, does contain enough quark-mass dependent struc- 
tures to reproduce the plateau-like behaviour at pion masses above 600 MeV without missing 
the physical value of gA- In other words, we suspect that the lattice data are just not ac- 
curate enough to properly constrain all three parameters g A , C SSE (1 GeV), g\ — which does 
not come as complete surprise given the low statistics of our data sample. In order to test 
this hypothesis we are trying to further constrain these couplings via available information 
from hadron phenomenology. An obvious idea is to fix the axial A A coupling g\ = 9/5 g A 
from SU(6) symmetry, leaving only g A and C as free parameters to be determined from the 
fit utilizing Eq. (fT2*|) . The numerical values resulting from this fit are shown in table El with 
the label fit II. The associated chiral extrapolation curve is shown in FigJ21 Its x 2 /d.o.f. is 
comparable with the values for the free fits la and lb. We conclude that — from a physical 
point of view — the SU(6) constrained fit does slightly better than the free fits when one 
considers the "higher" value for g A . Nevertheless, the axial coupling is still smaller than 1 
in the chiral limit, which is hard to reconcile with the Adler-Weisberger sum rule Eq.([Q. 
Surprisingly, even when an error analysis at the 95% confidence level is performed (see 
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Figure 2: The long-dashed curve represents the free fit to the QCDSF lattice data below 
m n = 750 MeV, utilizing the leading-one- loop SSE (Fit lb) results of Eq.(fT2"j). The analogous 
leading-one- loop HBChPT result (fit la) originating from Eq.fjHJ) is shown as the short- 
dashed curve. The solid line represents fit II for the SSE extrapolation with the additional 
SU(6) quark model constraint g\ = 9/5 g A . The indicated error band results from the 95% 
confidence ellipse shown in FigOl while the solid dot indicates the physical value of gA- 



FigHU), the resulting band shown in Fig. |2]does not include the physical point, despite the 
rather large error bars of the lattice data ! We thus conclude that the SU(6) assumption 
does not provide us with a reasonable mechanism that would allow for an enhancement of 
the chiral extrapolation curve around the physical pion mass. Similar conclusions can be 
drawn from the recent calculation by Detmold et al. [T^], which also utilized SU(6) assump- 
tions to reduce the number of unknown couplings. Interestingly, we note that the resulting 
value for C(lGeV) is still rather small, but has changed sign when compared to the free fit 
results. In the following we will pursue the hypothesis that it is the sign and magnitude 
of this coupling (and the resulting interplay with the non-analytic chiral results given in 
Eq.()12j0 that are responsible for the enhancement of gA{ m n) around physical pion mass. 

4.4 Fit Constrained by ttN — > ttttN 

When varying the three poorly known parameters g A , C(X), g\ of the SSE extrapolation 
formula Eq.([12|) one realizes that it is the short-distance physics encoded in C(A) which is 
crucial for a possible enhancement in gA^-n) at small pion masses. Empirically one finds 
that with negative values around for C(A = lGeV) ~ —4 GeV~ 2 it would be possible to 
counterbalance the trend of ^a(^tt) towards values below 1, which is driven by the chiral 
logarithm (cf. fit la in FigJ2J). The question arises whether such a value for C{\ = lGeV) 
at such a regularization scale is consistent with known physics. 
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Figure 3: Confidence ellipse for the SU(6)-constrained leading-one-loop SSE fit II shown in 
FigJ21 The two parameters shown denote the nucleon coupling constant g° A in the chiral 
limit versus the parameter C(A) discussed in the text, evaluated at the regularization scale 
A = 1 GeV. 



At first glance one would expect that only lattice simulations can provide more informa- 
tion on C(A) — respectively the couplings B 9 , B 2 o — to further constrain the fits. However, 
when analyzing the chiral tensors in Eq.(j3J), one notices that the physically observable cou- 
pling Bg also contributes in inelastic pion-nucleon scattering processes like ttN — > mrN 
|2~o] . In ref.j2H] this process was analyzed to 0(p 3 ) in HBChPT and values for the couplings 
of interest were obtained 5 . Here, however, we follow the substantially revised analysis of 
ref. |2H] • Translating these findings into our conventions we obtain 

B r 9 (\ = ml hys ) = (-1.4 ± 1.2) GeV" 2 , B r 20 (X = m^ 3 ) = , (14) 

with the error bar arising from three different fitting procedures described in j2H]- We note 
that refs. [2H1 12E] follow the convention to set all equation-of- motion dependent counterterms 
like B 2 q identically equal to zero at the scale of their analysis, which corresponds to A = 0.138 
GeV. While we can directly utilize this information in the heavy baryon analysis (c.f. the 
discussion in section EJ) there is no unique procedure how to import this information into the 
framework with explicit A(1230) degrees of freedom. The complication arises from the fact 
that Bg(X), -E?2o(A) have different beta-functions in HBChPT and SSE. One can therefore 
produce different chiral extrapolation curves depending on the scale where the matching 
between the two theories takes place. In the following we will perform the matching at the 
intermediate scale A = 2Ao = 0.54 GeV which suggests itself as a "natural" scale in SSE 
by looking at Eqs. ljlUlllJl . We have carefully examined that our results are not sensitive 
(within the uncertainties of the 7riV — > ttttN analysis) to changes of the scale A between 

5 A second analysis of irN — > ttttN scattering [22] unfortunately does not specify the value for the analogue 
of the coupling Bg they utilized. 
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0.4 and 0.8 GeV. At the end of this section we demonstrate that this matching choice is 
consistent with available information on the chiral limit behaviour of g^, which is the only 
constraint one can put on different matching prescriptions. First, we evaluate the couplings 
of Eq. (fT4|) at the scale A = 2 A via the relation 



i?[(2A ) = Bl(m?n - In , i = 9 , 20 (15) 



leading to the values 

B r g (X = 2A )\ HB = (-1.8 ± 1.2) GeV" 2 , B r w (X = 2A )\ HB = 0.9 GeV" 2 . (16) 



The leading-one-loop HBChPT beta-functions required can be found in ref. [20]. Note that 
we have utilized physical couplings in the heavy baryon beta-functions, inducing an uncer- 
tainty which is negligible compared to the error bar in Bg. At the scale A = 2A we then 
perform the matching 

B r 9 (\ = 2A )\ SSE = B r 9 (\ = 2A )\ HB , B r 20 (\ = 2A )\ SSE = B r 20 (\ = 2A )\ HB . (17) 

Next we run the SSE couplings up to the scale A = 1 GeV to facilitate an easier comparison 
with the previous fits (c.f. table |2J: 

3 SSE lGeV 

5[(lGeV)=5[(2A )--^--ln— — , z = 9, 20 . (18) 

y / iiTf 7T ) z 2A 



Utilizing the leading-one-loop SSE beta-functions of ref.[T4| with the couplings shown in fit 
III of table 121 one finally obtains 

B r 9 (X = lGeV)\ssE = (1.2 ± 1.2) GeV" 2 , B r 20 (X = lGeV)\ SS E = 1.9 GeV" 2 . (19) 

The error bar given for Sg(lGeV) originates from the determination of ref. [2H] shown in 
Eq. (jl4)l . No errors arising from the couplings in the HBChPT or SSE beta functions are 
shown, as they are smaller than the uncertainty in Bg resulting from the irN — ► inrN 
analyses. 

The resulting chiral extrapolation curve for gA{^n) utilizing Eqs. (jl2ll9|) and fitting 
g% gi to the five lattice data points is shown in FigEJ The resulting parameters are shown 
as fit III in table El The solid line in Fig0] refers to the central value for Bg in Eq. (jl9j) . 
whereas the two dashed curves indicate the influence of the (large) uncertainty about the 
precise value of Bg given in Eq. (J19|) . Nevertheless, while the x 2 /d.o.f. values are comparable 
with fits I and II shown in Fig|21 only fit III is consistent with the physical point g p J^ ys as 
shown in FigEJ We note that it is indeed the large negative values for the parameter 

C(lGeV) = 5 9 r (lGeV) -2g\ ££ (lGeV) = (-3.4 ± 1.2) GeV" 2 (20) 

that can produce an enhancement near the physical point. Given the large uncertainty in 
C(lGeV) (resulting from Bg) we believe that at present any normalization issues (e.g. result- 
ing from finite volume effects) can be accounted for by our extrapolation function ^a(^tt) of 
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Figure 4: Fit III: Incorporating known information from ttN — > ttttN differential cross 
sections into the leading-one-loop SSE analysis. The solid line shows the central value, 
whereas the dashed curves denote the upper and lower error bar from Eq.(|19|). The matching 
of the coupling constants was performed at A = 0.54 GeV as discussed in the text. 



Eq.()12|). Analyzing the ellipsoidal joint 95% confidence region for the fit parameters (g% g{) 
associated with the central value of C(lGeV) in fit III (c.f. Fig|SJ), the band of allowed 
fit curves completely covers the uncertainty region resulting from Bg (dashed curves) and 
shows that at this confidence level no fit solutions exist for which gA^n) becomes smaller 
than unity, in agreement with the Adler-Weisberger sum rule Eq.(Q). 

We therefore conclude that the strong downward-bending of gA^ir) driven by the Gold- 
stone boson cloud around spin 1/2 and spin 3/2 baryon states is counterbalanced by large 
short-distance corrections parameterized via C(A). The combination of these two effects 
leads to a rather flat chiral extrapolation function for gA{ m n), which — according to crucial 
input from ttN — > 7r7riV scattering analyses — can display an enhancement near the physical 
point to reconcile the lattice data with our knowledge from neutron beta decay. The analysis 
presented here suggests that in the (lattice) observable gA{ m Tr) no strong curvature is to be 
expected between the physical point and the location of the lattice data. The large LNA 
terms contained in gA{ m n) presumably only become visible for quark masses rather close to 
the chiral limit. 

Finally, we comment on the validity of our matching procedure resulting in the cou- 
plings of Eq.(Hnj). Once the couplings are specified at one particular scale the resulting 
SSE extrapolation curve shown in FigEJdoes not depend anymore on A, but it does depend 
on the matching point Eq. (jl7|) . Nevertheless, for each baryon observable there exists one 
unique curve through which the chiral limit value is approached when the quark masses go 
to zero. This is the reason why it is possible to utilize information about coupling constants 
in one version of ChEFT and transcribe it to another version of that theory. From Eq. fjll)) 
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Figure 5: Confidence ellipses for the constrained leading-one-loop SSE fit III shown in FigHJ 
The 2 parameters shown denote the nucleon coupling constant g A in the chiral limit versus 
the axial A A coupling g\. The three ellipses shown correspond to central, upper and lower 
values for the short distance couplings determined in Eq. (fT9|) at a scale of A = 1 GeV. 



we already know that the first two terms in the chiral expansion of 5u(^tt) are identical in 
HBChPT and SSE, providing an important consistency check. However, the terms ~ can 
be different between the two ChEFTs — even when each one is evaluated at leading-one-loop 
order as done here — because the two schemes can build up the contributions to a particular 
term in the chiral expansion more or less effectively in their respective perturbative series. 
Utilizing Eq. (jl6|) one finds 

77? 

g\\ HB = l.2 + 4 1 HB ml\n-^ + ml(-l7.8±4.8)GeV- 2 + 0(ml) (21) 

for the chiral limit behaviour of the leading-one- loop HBChPT result of Eq.fjHJ). On the 
other hand, utilizing the parameters of fit III (c.f. table EJ) one obtains 

77? 

g A \ SSE = 1.2 + 4 7 HB m>-f- + (-10.4 ± 4.8)GeV^ 2 + O^) (22) 

for the chiral limit behaviour of the leading-one-loop SSE result of Eq. (fT2*|) . The m\ terms 
agree within the error bars and the shifts of the central values could reasonably be accounted 
for by higher order corrections to the heavy baryon result of Eq. (J21|) . We conclude that the 
matching procedure used in this work is consistent with all presently known information 
about the chiral limit of su^ti-)- In the next section we will test the stability of fit III by 
including the physical point as an additional constraint for the fit and then conclude our 
analysis. 
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5 Discussion 



In the previous section we have found that the leading-one-loop SSE result Eq.([T2|) together 
with input from scattering analyses can provide a meaningful chiral extrapolation of the 
present lattice data which is consistent with the experimentally known value for g^. In 
order to test the stability of our result we now study the changes to our curves when we 
include the information at the physical point as an additional constraint together with the 
lattice data. Reassuringly, one finds only small changes in the SSE case. Numerical values 
for the set of couplings labeled fit IVb are given in table El entirely consistent with the 
couplings of fit III which does not know about the location of the physical point. The 
resulting extrapolation function is shown in FiglHJ providing a sensible connection between 
the chiral limit value, the physical point and the QCDSF lattice data. We note that a 
large negative value of C(lGeV) and g% > gi U ^ are common features of fits III and IVb, 
indicating that our reasoning provided in the previous section was on the right track. 

When taking the physical point as additional input, we can now also examine the validity 
of the leading-one-loop HBChPT extrapolation formula of Eq.fjSJ). In fit IVa we have taken 
C(lGeV) = —5.3 GeV -2 , which corresponds to the (central) value suggested in ref . . 
The resulting (short-dashed) curve shown in FigJU] completely misses the lattice data and 
leads to a substantially different chiral limit value g% even if the rather large uncertainty in 
C(lGeV) (c.f. table|21) — indicated via the bands in Fig|7] — is taken into account. According 
to our analysis this finding forces us to the conclusion that the leading-one-loop HBChPT 
curve as given in Eq.fjHJ) is not sufficient to even perform the chiral extrapolation from the 
chiral limit to the physical point! With the help of the decoupling constraints Eq. ljlO)) we 
have made sure that the chiral limit value of ^(^tt) is the same in both approaches, as there 
exists only one definite axial charge of the nucleon in the chiral limit. In order to further test 
the applicability of the HBChPT curve we have also performed fit IVc, demanding that g° A 
takes the same value as in the leading-one-loop SSE analysis. The resulting (long-dashed) 
extrapolation curve is shown in FiglHl It only follows the SSE curve for extremely small 
values of m„ misses the physical point and breaks down before it comes anywhere close to 
the existing lattice data. We therefore consider the use of Eq.(jHJ) as not appropriate for a 
stable chiral extrapolation outside the region of the extreme chiral limit. 

We note that it is the very leading-one- loop HBChPT formula Eq.fjHJ) that has been 
utilized recently in refs. |29l I3*U] to study the behaviour of the deuteron in the chiral limit. 
In view of the findings presented here the question arises whether conclusions drawn in 
refs. [23 HID] regarding the bound state nature of the deuteron are unaltered, when the chiral 
extrapolation of ^("v) from the physical point down to the chiral limit follows the curve 
from fit IVb, rather than the one from IVa as shown in FigEl 

6 Summary and Outlook 

To conclude, we summarize the main results of our work: 

- The leading-one-loop HBChPT formula for ^(m^) Eq.(JE} depends on two free pa- 
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Figure 6: Fit IV: Incorporating the physical point and the lattice data into the leading-one- 
loop SSE extrapolation formula results in a well behaved chiral extrapolation function (solid 
line, fit IVb) from the chiral limit over the physical point out to the region of the lattice 
data. The short-dashed line shows the corresponding leading-one-loop HBChPT result (fit 
IVa). The long-dashed line (fit IVc) indicates the HBChPT result assuming the same chiral 
limit value g A as in the SSE analysis of fit IVb. 



rameters g A , C(A). Crucial information on sign and size of C(A) can be obtained from 
analyses of the nN — > ttttN reaction. C(A) is large and negative for typical values of 
the regularization scale. 

- The leading-one- loop SSE formula for gA{mn) Eq. (|12jl can be matched to the HBChPT 
formula in the chiral limit. It has sufficient quark mass dependent structures to pro- 
vide a consistent chiral extrapolation from the chiral limit across the physical point up 
to the lattice data, assuming that the large quark mass quenched lattice data are con- 
sistent with fully dynamical simulations and that finite volume effects do not provide 
large corrections. 

- Present lattice data are not precise enough in order to determine all couplings of the 
SSE formula from the data alone. Additional information from ttN dynamics has to be 
invoked. SU(6) arguments are not sufficient to generate the required form of gjijjn^). 

- The chiral extrapolation curve for gA{ m n) is rather flat. It implies a co-operation of 
terms in the quark mass expansion such that pion cloud effects are basically "turned 
off" for pion masses exceeding 0.6 GeV. A small enhancement is seen around the phys- 
ical pion mass. This enhancement results from an interplay between the Goldstone 
boson dynamics of the nucleon and short distance contributions which can be con- 
strained from ttN — > ttttN scattering data. The chiral limit value of the axial coupling 
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Figure 7: We have magnified the region around the chiral limit up to the physical point. 
The bands indicated for the two heavy baryon curves arise from the uncertainty of the 
Ti N — > Tin N analysis of Eq.(JTH). 



of the nucleon is found to be 



^« 1.2 ±0.1. (23) 

The parameters of fit IVb give the best available description of our present knowledge 
on gA{m n ). 

- The leading one-loop HBChPT formula (JHJ) for gA{™-x) is not useful for chiral extrap- 
olation. It fails to connect the chiral limit with the physical value of the pion mass. 
According to our analysis the chiral logarithm as the leading non-analytic quark mass 
term in the HBChPT formula is only visible in extrapolations performed below the 
physical point. 

Finally, we note that the role of finite size effects in the lattice data for g& has to be analyzed 
in a systematic framework to better understand the systematics of the lattice error bars. 
More data from simulations for ^(m^) performed at different lattice volumes are called for 
to guide these studies. We are confident that ChEFT is an important tool to resolve this 
issue in the future. 
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A Amplitudes 



Here we present the results for the 8 leading-one-loop amplitudes, i.e. 0(e 3 ) diagrams 
shown in FigC] which can contribute 6 to the axial coupling constant g A of a nucleon. The 
Lagrangians needed for this calculation are discussed in the main text. We work in the 
Breit-frame and choose the velocity vector v^ 1 = (1,0,0,0). With S* M denoting the Pauli- 

-A 

finds 



Lubanski spin- vector and e A the polarization 4- vector of an external axial- vector source, one 
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(24) 

(25) 
(26) 
(27) 

(28) 



Explicit expressions for the function J 2 are given in [31] . We evaluate the amplitudes in d 
dimensions with induced regularization scale A. Any ultraviolet divergences appearing in 
the limit d — > 4 are subsumed in 



A'' 



16tt 2 



^ + 1(7,-1 -Into) 



(29) 



where 7^ denotes the Euler-Mascheroni constant. 

Finally, we note that with the decoupling choice made in Eq. (JlOJ) the nucleon Z-factor 
calculated to leading-one-loop in SSE reads 



7SSE 
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G(e 4 



(30) 



where i? has been defined in Eq. (ll3|) . 



6 To simplify the calculation we have made use of the gauge-condition e A • q = 0. 
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